clear;
clc;
A=[7,5,5,8;
    5,6,9,7;
    5,9,9,0;
    8,7,0,1];
[D,p]=chol(A)%代码中的p用来判断该矩阵是否是正定矩阵